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We perform a systematic study on the surface property of nucleus-nucleus potential in heavy-ion 
reactions using large-angle quasielastic scattering at energies well below the Coulomb barrier. At 
these energies, the quasielastic scattering can be well described by a single-channel potential model. 
Exploiting this fact, we point out that systems which involve spherical nuclei require the diffuseness 
parameter of around 0.60 fm in order to fit the experimental data, while systems with a deformed 
target between 0.8 fm and 1.1 fm. 
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I. INTRODUCTION 

The Woods-Saxon form, which is characterized by 
the depth, radius and diffuseness parameters, has often 
been used for the inter-nuclear potential for heavy-ion 
reactions. Conventionally, the diffuseness parameter of 
around 0.63 fm has been employed for calculations of 
elastic and inelastic scattering, which are sensitive only 
to the surface region of the nuclear potential P, P. This 
value of surface diffuseness parameter has been well ac¬ 
cepted, partly because it is consistent with a double fold¬ 
ing potential 0. In contrast, a recent systematic study 
has shown that experimental data for heavy-ion fusion 
reactions at energies close to the Coulomb barrier re¬ 
quire a larger value of the diffuseness parameter, ranging 
between 0.75 and 1.5 fm, as long as the Woods-Saxon 
parameterization is used as a nuclear potential 0. The 
origin of the discrepancy in the surface diffuseness pa¬ 
rameter between the scattering and fusion processes has 
not yet been understood. 

Large-angle quasielastic scattering at deep sub-barrier 
energies provides an alternative way to look at this prob¬ 
lem. Quasielastic scattering and fusion are both in¬ 
clusive processes and are complimentary to each other. 
The former is related to the reflection probability at 
the Coulomb barrier, while the latter to the penetra¬ 
tion probability. In heavy-ion reactions at energies near 
the Coulomb barrier, it is well known that the channel 
coupling effects caused by the collective inelastic excita¬ 
tions of the colliding nuclei strongly affect the reaction 
dynamics At deep sub-barrier energies, however, 

the channel coupling effects on quasielastic scattering can 
be disregarded, since the reflection probability is almost 
unity at these energies irrespective of the presence of 
channel couplings, even though inelastic channels them¬ 
selves may be strongly populated |3]- This is similar to 
fusion at energies well above the Coulomb barrier, where 
the penetrability is almost unity 0. 

The above concept was recently applied to the exper¬ 
imentally measured quasielastic scattering cross sections 


for the -I- ^^^Sm system at deep sub-barrier energies 
0. It was found that the larger surface diffuseness pa¬ 
rameter of around 1.0 fm is required for this system in 
order to fit the data. This value is consistent with the 
one required for fusion. 

It is apparent that a more systematic study is nec¬ 
essary, in order to clarify whether the quasielastic scat¬ 
tering around the Coulomb barrier generally requires a 
larger value of surface diffuseness parameter than the 
conventional value of around 0.63 fm. The aim of this 
paper is to carry out such systematic study on quasielas¬ 
tic scattering at deep sub-barrier energies. To this end, 
we calculate the excitation function of the quasielas¬ 
tic cross sections for systems involving both spheri¬ 
cal and deformed nuclei. The reactions ^^’^'*S-|-^®^Au, 

32.34g_^208p,3^ 16Q_^154gj^^ 186y^^ 208pp^ foj. g^- 

perimental data exist at deep sub-barrier energies, are 
studied. We show that a surface diffuseness parameter 
of around 0.6 fm is favored by the data for reactions in¬ 
volving spherical nuclei, whilst those involving deformed 
nuclei require a larger value of the diffuseness parameter. 

The paper is organized as follows. In the next section, 
we briefly review the large-angle quasielastic scattering 
at deep sub-barrier energies. We also explain the pro¬ 
cedure of our analyses which use a one dimensional ion- 
ion potential, including our definition of deep sub-barrier 
energies. In Sec. Ill, we present our results for the 
fitting and discuss its sensitivity to the barrier height en¬ 
ergy and to the channel coupling effects. We summarize 
the paper in Sec. IV. 


II. METHOD OF ANALYSES 

A. Large-angle quasielastic scattering at deep 
sub-barrier energies 

Our purpose in this paper is to study the surface prop¬ 
erty of ion-ion potential using heavy-ion quasielastic scat¬ 
tering. Before we explain the method of our analyses, let 
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us first discuss briefly the advantage of exploiting large- 
angle quasielastic scattering at deep sub-barrier energies. 

At energies well below the Coulomb barrier, the cross 
sections of (quasi)elastic scattering are close to the 
Rutherford cross sections, with small deviations caused 
by the effect of nuclear interaction. This effect can be 
taken into account by the semiclassical perturbation the¬ 
ory. The ratio of elastic scattering Uei to Rutherford cross 
sections an at a backward angle 9 is given by H, 

daeijE 

cm 1 ^ 1 _l_ Cjv(re) V2a7rfc7? 

9') kcL Eqjyi 

where Ecm is the centre-of-mass energy, k = ^/2jlE^/h, 
/r being the reduced mass, and rj is the Sommerfeld pa¬ 
rameter. This formula is obtained by assuming that 
the nuclear potential Vv('c) has an exponential form, 
exp(—r/a), around the classical turning point rc = {ij + 
+ A 2 )/fc, where Ac = r?cot( 0 / 2 ) is the classical an¬ 
gular momentum for the Rutherford scattering. We see 
from this formula that the deviation of the elastic cross 
sections from the Rutherford ones is sensitive to the sur¬ 
face region of the nuclear potential, especially to the sur¬ 
face diffuseness parameter a. Notice that, for small scat¬ 
tering angles, the Fresnel oscillation may complicate the 
formula. Also, as mentioned in the previous section, the 
channel coupling effects on the quasielastic cross sections 
are negligible at deep sub-barrier energies. We can thus 
study the effect of the surface diffuseness parameter in a 
transparent and unambiguous way using the large-angle 
quasielastic scattering at deep sub-barrier energies. 


B. Procedure 

In order to compare with the experimental data for 
the quasielastic cross sections at deep sub-barrier ener¬ 
gies, we use a one-dimensional optical potential with the 
Woods-Saxon form. Absorption following transmission 
through the barrier is simulated by an imaginary poten¬ 
tial with W = 30 MeV, = 0.4 fm, and = 1.0 fm. 
This model calculates the elastic and fusion cross sec¬ 
tions, in which the elastic cross sections can be considered 
as quasielastic cross sections to a g ood approximation at 
these deep sub-barrier energies [13 . Note that the results 
are insensitive to the parameters of the imaginary part 
as long as it is well localized inside the Coulomb barrier. 

In order to carry out a systematic study, we calculate 
the Coulomb barrier height using the Akyiiz-Winther po¬ 
tential M- We examine several potentials with different 
values of surface diffuseness parameter, which give the 
same calculated barrier height. To this end, we vary the 
radius parameter rg while keeping the depth parameter 
Vo to be 100 MeV. This is possible because the effect 
of variation in Vq and rg on the Coulomb barrier height 
compensates with each other at the surface region. 

We define the region of “deep sub-barrier energies” in 
the following way. In heavy-ion collisions at energies near 


the Coulomb barrier, collective inelastic excitations of the 
colliding nuclei and transfer reactions are strongly cou¬ 
pled to the relative motion. This causes the splitting 
of the Coulomb barrier into several distributed barriers 
mil. We define the deep sub-barrier energies as around 
3 MeV below the lowest barrier height or smaller. For this 
purpose, we first use the computer code CCFULL ^ in 
order to explicitly construct the coupling matrix (which 
includes the excitation energy for the diagonal compo¬ 
nents) for the coupled-channels equations for each sys¬ 
tem by including known low-lying collective excitations. 
We then diagonalize it to obtain the lowest eigen-barrier. 

We find that the deep sub-barrier region defined in 
this way corresponds to the region where the experimen¬ 
tal value of the ratio of the quasielastic to the Ruther¬ 
ford cross sections is larger than around 0.94. We there¬ 
fore include only those experimental data which satisfy 
d<7qe\/duR > 0.94 in the fitting. A few experimental 
data points with values exceeding unity were excluded 
while performing the fits, but are shown in the figures 
below. 

We apply this procedure to the 32,34g _|_ 208 p| 2 ^ 32,34g 
-I- ^®^Au 01, and -I- ^°®Pb ^3 reactions which in¬ 
volve spherical nuclei, as well as the -I- ^^^Sm and 
-I- ^®®W reactions ^3 which involve a deformed tar¬ 
get. For the deformed systems the scarcity of data points 
at deep sub-barrier energies led us to extend the fitting 
region to somewhat higher energies. This meant that the 
calculations had to take account of deformation effects as 
explained in Sec. III. B. 


III. RESULTS AND DISCUSSION 
A. Spherical systems 

We first present the results for systems involving spher¬ 
ical nuclei. Figure 1 compares the experimental data 
with the calculated cross sections obtained with different 
values of the surface diffuseness parameter in the Woods- 
Saxon potential for the -I- ^®’^Au system (the upper 
panel) and the -I- ^®^Au system (the lower panel). 
The Coulomb barrier height is 141.2 MeV for the -I- 
^®^Au reaction and is 140.2 MeV for the -I- ^®^Au re¬ 
action. The best fitted values for the surface diffuseness 
parameter are a = 0.57 ± 0.04 fm and a = 0.53 ± 0.03 fm 
for the and -I- ^®'^Au reactions, respectively. The 
cross sections obtained with these surface diffuseness pa¬ 
rameters are denoted by the solid line in the figure. The 
dotted and the dot-dashed lines are calculated with the 
diffuseness parameter of a = 0.80 fm and a = 1.00 fm, re¬ 
spectively. Figure 2 shows the results for the -I- ^°®Pb 
(the upper panel) and the -I- ^°®Pb (the lower panel) 
reactions. The Coulomb barrier height is 145.1 MeV and 
144.1 MeV for the and -I- ^°®Pb reactions, respec¬ 
tively. The best fitted values for the surface diffuseness 
parameter are a = 0.60 ± 0.04 fm and a = 0.63 ± 0.04 fm 
for the and -I- ^°®Pb reactions, respectively. 
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FIG. 1: The ratio of the quasielastic to the Rutherford cross 
sections at Ouh = 159° for the + ^®^Au (the upper panel) 
reaction and for the + ^®^Au (the lower panel) reaction. 
The experimental data are taken from Ref. Hi-' The solid 
line results from using a diffuseness parameter obtained by 
performing a least-square fit to the data. The dotted and the 
dot-dashed lines are obtained with the diffuseness parameter 
of a = 0.80 fm and a = 1.00 fm, respectively. 


It is evident from Figs. 1 and 2 that these spherical 
systems favor the standard value of the surface diffuse- 
ness parameter, around a = 0.60 fm. The calculations 
with the larger diffuseness parameters, a = 0.80 fm and 
1.00 fm, underestimate the quasielastic cross sections and 
are not consistent with the energy dependence of the ex¬ 
perimental data. We obtain a similar conclusion for the 
isq _|_ 208 p|^ system, where the best fitted value for the 
surface diffuseness parameter is a = 0.59 ± 0.10 fm with 
the Coulomb barrier height of 76.1 MeV. The result for 
this system is shown in Fig. 13 

The conclusions are not sensitively dependent on the 
choice of barrier height energy Vb ■ In order to demon¬ 
strate this, we vary the barrier height by 1%, and repeat 
the same analyses. The result for the -I- ^°®Pb system 
is shown in Fig. 2| The solid line denotes the result ob¬ 
tained with the Akyiiz-Winther potential, as a reference, 
which is the same as the solid line in the upper panel of 
Fig. 2. The best fits and the resulting a values using Vb 
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FIG. 2: The ratio of the quasielastic to the Rutherford cross 
sections for the -I- ^°®Pb (the upper panel) reaction at 
^lab = 170° and for the ^“^8 -I- ^°®Pb (the lower panel) reaction 
at Slab = 159°. The experimental data are taken from Ref. 
H. The meaning of each line is the same as in Fig. 1. 
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FIG. 3: The ratio of the quasielastic to the Rutherford cross 
sections for the ^®0 -I- ^®®Pb reaction at Slab = 170°. The 
experimental data are taken from Ref. m. The meaning of 
each line is the same as in Fig. 1. 
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FIG. 4: Comparison of quasielastic cross sections obtained 
for three different values of the Coulomb barrier height for 
the + ^°®Pb reaction. The surface diffuseness parameter is 
determined for each barrier energy by htting the experimental 
data. 
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FIG. 5: The ratio of the quasielastic to the Rutherford cross 
sections for the + ^®^Sm reaction at 0iab = 170°. The 
solid line is obtained using the best htted value of the surface 
diffuseness parameter, a = 1.14 fm. The dotted line denotes 
the cross sections obtained with the diffuseness parameter of 
a = 0.60 fm. The experimental data are taken from Ref. M- 


= 143.6 MeV and Vb = 146.5 MeV are also shown in Fig. 
21 The a value changes by ±0.04 fm for a ±1% change 
in the barrier energy. The cross sections obtained with 
these potentials are shown in the hgure by the dotted 
and the dot-dashed lines, respectively. One clearly sees 
that the effect of the variation of the Coulomb barrier 
height on the surface diffuseness parameter is small. The 
barrier energy obtained from the analysis of the above¬ 
barrier fusion cross sections is 144.03 MeV 0, which is 
within the range of Vb used in the calculations. Thus, 
the diffuseness parameter extracted in this work will not 
change signihcantly if Vb determined from fusion data, 
instead of the Akyiiz-Winther prescription, is used. We 
have conhrmed a similar behavior of the surface diffuse¬ 
ness parameter a for the other systems as well. 

B. Deformed systems 

Let us next discuss the systems with a deformed tar¬ 
get, that is, ± reactions. For these 

systems, only a few data points are available at deep sub¬ 
barrier energies. We therefore include the experimental 
data at energies not only well below but also around the 
lowest barrier in the htting procedure. At these en¬ 
ergies, the channel coupling effects start playing an im¬ 
portant role in quasielastic reactions, and we include the 
effect of deformation of the target nucleus in our calcula¬ 
tions. Therefore, our analyses for the deformed systems 
are somewhat more model dependent than those for the 
spherical systems presented in the previous subsection. 

In order to account for the deformation effect on the 
quasielastic scattering, we use the orientation average for¬ 
mula in which we neglect the Hnite excitation 


energy of the ground state rotational band. With this 
formula, the quasielastic cross section is given by. 


^qel (A'cm ; — / d(cOS (Tel (F^cm j ^t)5 (2) 

Jo 

where 9t is the angle between the symmetry axis of the 
deformed target and the direction of the projectile from 
the target. In the calculation for both the systems, we 
take six different orientation angles into account [l^ . 
The results change only marginally even if we include 
the larger number of orientation angles. 

The best fitted value for the surface diffuseness param¬ 
eter obtained in this way is a = 1.14 ± 0.03 fm and 0.79 
±0.04 fm for the ± and ± ^®®W reactions, 

respectively. The deformation parameters which we use 
in the calculations are (32 = 0.306 and P 4 = 0.05 for 
^^"‘Sm and P 2 = 0.29 and P 4 = —0.03 for ^®®W. Figs. El 
and El compare the calculated cross sections with the ex¬ 
perimental data. The solid line in each hgure is obtained 
using the best htted value of the dihuseness parameter. 
The dotted line shows the cross section obtained with the 
dihuseness parameter of o = 0.60 fm as a reference. We 
hnd that the larger values of the surface dihuseness pa¬ 
rameter, a = 1.14 fm and 0.79 fm, in the nuclear potential 
are favored for these deformed system, in accordance with 
our previous conclusion in Ref. Q For the ± 
reaction, the calculated cross sections with the standard 
value of the surface dihuseness parameter around 0.60 fm 
are clearly in disagreement with the experimental data. 
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FIG. 6: The ratio of the quasielastic to the Rutherford cross 
sections for the reaction at 0iab = 170°. The 

meaning of each line is the same as in Fig. |3 The experi¬ 
mental data are taken from Ref. IT^ . 


C. Discussion 
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FIG. 7: The best fitted values of the surface diffuseness pa¬ 
rameter a as a function of the charge product of the projectile 
and target nuclei, ZpZt. The filled circles and triangles are 
for the spherical and the deformed systems, respectively. The 
open circles and triangles are the surface diffuseness parame¬ 
ters deduced from the analyses of fusion cross sections Q. 



Figure [ 7 | summarizes the results for our systematic 
study for the surface diffuseness parameter. It shows the 
best fitted value of diffuseness parameter as a function of 
the charge product of the projectile and target nuclei for 
each system. The results for the spherical systems are de¬ 
noted by the filled circles, while those for the deformed 
systems the filled triangles. One clearly sees the trend 
that the best fitted value of the diffuseness parameter is 
around 0.60 fm for the former, while it is much larger 
than that for the latter. Also, one sees that the surface 
diffuseness is almost constant for the spherical systems. 

The value of the surface diffuseness parameter obtained 
in this study for the spherical systems agrees well with 
the conventionally used value a ~ 0.63 fm. This suggests 
that the double folding potential is valid at least in the 
surface region and for systems which do not involve a 
deformed target. For these systems, the discrepancy be¬ 
tween the values of the diffuseness parameter determined 
from fusion data (open circles and triangles in Fig. [71) 
and those from quasielastic data must be related with 
the dynamics inside the Coulomb barrier 0. 

For the deformed systems studied here, the diffuse¬ 
ness parameter extracted from the quasielastic scatter¬ 
ing is much larger than the conventional value of a ^ 
0.63 fm. Although this value is consistent with that ex¬ 
tracted from fusion, the origin of the difference between 
the spherical and the deformed systems is not clear. One 
should bear in mind, however, that our analyses for the 
deformed systems are somewhat model dependent. This 
is due to the fact that the experimental data in the deep 
sub-barrier region are sparse for the deformed systems, 
and we need to include the deformation effect in the cal¬ 
culations in order to reproduce the strong energy depen¬ 


dence of the quasielastic cross sections at energies around 
the lowest barrier where the data exist. In order to clar¬ 
ify the difference in the diffuseness parameter between 
the spherical and the deformed systems, further preci¬ 
sion measurements for large-angle quasielastic scattering 
at deep sub-barrier energies will be necessary, especially 
for deformed systems. 


IV. SUMMARY 

Large-angle quasielastic scattering provides a powerful 
tool not only for the analysis of the barrier distribution 
around the Coulomb barrier but also for the study of the 
surface property of the nuclear potential. This is due to 
the fact that channel coupling effects play a minor role in 
quasielastic scattering at deep sub-barrier energies, that 
enables a relatively model independent analysis of ion- 
ion potential. Using this fact, we have systematically 
analyzed experimental data for quasielastic scattering at 
deep sub-barrier energies, with the aim of extracting the 
surface diffuseness parameter of internuclear potential. 
We obtained the diffuseness parameter that is consistent 
with the standard value of around a = 0.63 fm for the 
systems involving spherical nuclei. In contrast, fits to 
the data for systems involving deformed nuclei require 
diffuseness parameter to be in the range of 0.8 to 1.1 fm, 
similar to that obtained from analyses of the fusion data 
at above-barrier energies. 

The origin of the difference between the spherical and 
the deformed systems is not clear at the moment. In 
order to clarify this and confirm the systematics found 
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in this paper, more experimental investigations on large- 
angle quasielastic scattering at deep sub-barrier energies 
will be certainly helpful, especially for deformed targets. 
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